Semiclassical theory of the Ehrenfest-time dependence of quantum transport 
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In ballistic conductors, there is a low-time threshold for the appearance of quantum effects in 
transport coefficients. This low-time threshold is the Ehrenfest time te. Most previous studies of 
the TE-dependence of quantum transport assumed ergodic electron dynamics, so that they could 
be applied to ballistic quantum dots only. In this article we present a theory of the TE-dependence 
of three signatures of quantum transport — the Fano factor for the shot noise power, the weak 
localization correction to the conductance, and the conductance fluctuations — for arbitrary ballistic 
conductors. 
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I. INTRODUCTION 

The fact that electrons are described by quantum me- 
chanics, not classical mechanics, can be confirmed from 
a number of features in the electrical transport proper- 
ties of metals at low temperatures Well-known exam- 
ples of such 'signatures of quantum transport' are the 
Fano factor for the shot noise power, the weak localiza- 
tion correction to the conductance, and the conductance 
fluctuations of a mesoscopic conductor. Disorder plays 
an important role in the understanding of these effects 
in normal metals, because electrical transport in metals 
at low temperatures is dominated by scattering off im- 
purities comparable in size to the Fermi wavelength of 
the conduction electrons. Theoretically, the presence of 
disorder allows for the use of powerful field-theoretic tech- 
niques, which give accurate predictions for an 'ensemble 
average' over different impurity configurations.^'^ 

An important feature of the signatures of quantum 
transport in disordered conductors is that their size is in- 
dependent of material properties, sample size, or the con- 
centration or type of impurities lii^ii^ii^i^ The conditions 
under which this 'universality' appears are rather mild: 
At zero temperature and for dc transport one needs^ 



(1) 



where Ai? is the Fermi wavelength, vf the Fermi veloc- 
ity, r the elastic scattering time, and td the typical dwell 
time of electrons travelling between source and drain con- 
tacts. In addition one requires that the sample's dimen- 
sionless conductance g ^ 1, so that Anderson localiza- 
tion can be ruled out. The quantum effects do, however, 
show a weak dependence on the nature of the electron dy- 
namics near td, which is why the 'universal' signatures 
of quantum transport have slightly different magnitudes 
in, e.g., disordered quantum dots and quantum wires, re- 
flecting the difference between ergodic and diffusive dy- 
namics in these two systems. 

With the fabrication of high-mobility two-dimensional 
electron gases in semiconductor heterostructes, it has be- 
come possible to study quantum transport in devices in 



which the electron motion is ballistic over significantly 
longer distances than in metal samples, without scatter- 
ing off point-like impurities. In such devices, nontriv- 
ial geometries are achieved by the placement of artificial 
scattering centers or sample boundaries. Two paradig- 
matic examples, a ballistic quantum dot and an antidot 
lattice, are shown in Fig.[T] A quantum dot is a region of 
a two-dimensional electron gas confined by metal gates 
and coupled to source and drain electrodes via narrow 
contacts;^ An antidot lattice is an electron gas with arti- 
ficial macroscopic scattering centers.-* In the theoretical 
literature, a collection of randomly placed circular anti- 
dots is referred to as a 'Lorentz gas'. 

Because the signatures of quantum transport in disor- 
dered metals do not depend on impurity concentration or 
type one may be tempted to expect that ballistic conduc- 
tors are characterized by the same universal signatures of 
quantum transport as their disordered counterparts. The 
equivalent expectation in the context of spectral statistics 
is known as the "Bohigas-Giannoni-Schmit conjecture"— 
and believed to be true. In a seminal article, Aleiner 
and Larkin pointed out that this expectation need not 
always be correct for transport, however. They argued 
that there is a minimal time required for the appearance 
of quantum effects in ballistic conductorSi ^^i^^ This time 
is the 'Ehrenfest time' te, the time it takes for a minimal 
wavepacket to diverge and reach a size such that it can no 
longer be described by a single classical trajectory i^°i^^ 
The Ehrenfest time poses a short-time threshold for the 
appearance of quantum effects, because quantum phe- 
nomena cannot occur as long as a wavepacket travels 
along a single classical trajectory. One expects the same 
signatures of quantum transport in disordered and bal- 
listic conductors only if te ^ td-^'^ 

For a ballistic conductor in which the classical electron 
dynamics is chaotic with Lyapunov exponent A, one has 



TE = ^ \n{Ls/XF), 



(2) 



where Lg is a classical separation beyond which trajec- 
tories should be considered uncorrelated. In most ex- 
periments, A ~ and the logarithm in Eq. ^ not 
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numerically large, so that te is not much larger than 
the elastic scattering time r. This explains why, in- 
deed, the signatures of quantum transport in many bal- 
listic conductors are so similar to those of normal met- 
als with point-like impurities. ■^■^^ Nevertheless, there is 
no a priori reason why the logarithm in Eq. ([2]) must 
be small, and one may ask about the fate of quan- 
tum transport if te ~ td (or te comparable to in- 
verse frequency or the appropriate inelastic time, if time- 
dependent transport or finite temperatures are consid- 
ered). This question has received increasing attention in 
the last decade^i^^'ii^'ia^'Si^^^^^^^^^ 
It is a question of fundamental importance from a theo- 
retical point of view, because the regime of large te is the 
only parameter regime in which the signatures of quan- 
tum transport may discriminate between disordered and 
ballistic conductors. Moreover, large Ehrenfest times ap- 
pear naturally in the semiclassical limit Xp vpT, which 
provides one of the conditions necessary for the univer- 
sality of the signatures of quantum transport, see Eq. ([T]) 
above. 

The Ehrenfest-time dependence of weak localization 
was first addressed in the original article by Aleiner and 
LarkinpiS. The theory of Ref. [l^ is based on a field- 
theoretic approach which employs a minimal amount of 
diffraction from disorder in order to mimic the diffractive 
effects of scattering off the curved boundaries of the sam- 
ple and the artificial scattering centersi^ For the Lorentz 
gas, Aleiner and Larkin showed that the weak localization 
correction 5a to the ac conductivity acquires Ehrenfest- 
time dependent oscillations, 5a oc exp(2iLi;rE). They also 
considered the weak localization correction 5G to the dc 
conductance of a ballistic quantum dot, which is pro- 
portional to exp(— te/td — TE/T0)f2^ being the dot's 
dephasing time. Agam, Aleiner, and Larkin calculated 
the Fano factor / of a ballistic quantum dot,^^ which has 
the same exponential dependence cx exp(— te/td) as 5G. 
The experimental observation of the suppression of weak 
localization in an antidot lattice at large t-^/t^ and the 
suppression of the shot noise power in a ballistic quantum 
dot at large t^/to were consistent with the theoretical 
predictions fi^iii 

There is a second theoretical approach to quantum 
transport in ballistic conductors. This approach starts 
from a semiclassical expression of the sample's scattering 
matrix in terms of classical trajectories connecting the 
contacts. Since the conductance is proportional to the 
square of a scattering amplitude, the conductance is then 
expressed as a double sum over classical trajectories.'^^ 
Originally, the trajectory sums were performed in the so- 
called diagonal approximation, in which only the diago- 
nal terms in the double sum were kept.^i^i^l Although 
the diagonal approximation could explain the existence of 
weak localization and universal conductance fluctuations 
in ballistic quantum dots, as well as the dependences on 
the Fermi energy and an applied magnetic field. 
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could not describe the Ehrenfest-time dependences. 
A technical breakthrough occurred when Sieber and 



FIG. 1: Schematic picture of a ballistic quantum dot (left) 
and a Lorentz gas (right). In a ballistic quantum dot, the 
conductance is dominated by the conductances of the two 
contacts that connect the quantum dot to the electron reser- 
voirs. A Lorentz gas is a ballistic conductor with scattering 
from circular discs. 



Richter were able to include the leading off-diagonal 
terms into the summation. '^^'■^^ With off-diagonal con- 
tributions, the trajectory-based approach could success- 
fully capture te dependences, not only of the weak 
localization correction and Fano factor of a quantum 
dot j-'^^i^^'^'^i^^'^^i^^ but also of quantum signatures whose 
Te dependence was not known from the field-theoretic 
approach, such as the conductance fluctuations or the 
current pumped through a quantum dot with time- 
dependent shape (a "quantum pump")i^i^i^ Unlike 
weak localization and the Fano factor, the variance of the 
conductance var G and the mean square pumped current 
were found not to disappear in the limit of large Ehrenfest 
times. In fact, in the absence of dephasing, var G is inde- 
pendent of Te in a quantum dot, as was first observed by 
Tworzydlo et al. and Jacquod and Sukhorukov on the ba- 
sis of numerical simulations.^^ In the semiclassical theory 
the remarkable TE-insensivity of var G in quantum dots 
has its origin in a large contribution to var G from trajec- 
tories that spend a long time in the vicinity of periodic 
orbits. Since such trajectories can have arbitrarily long 
dwell times, the existence of te as a short-time threshold 
for quantum effects no longer poses a limitation on the 
size of mesoscopic fluctuations. 

With the exception of the original article by Aleiner 
and Larkin,!- who considered weak localization in a 
Lorentz gas, all theoretical work on te dependences has 
focused on ballistic quantum dots. The goal of the 
present article is to investigate the te dependence of 
quantum transport in arbitrary ballistic conductors. We 
use the trajectory-based approach and restrict ourselves 
to dc transport at temperatures low enough that dephas- 
ing does not play a significant role, so that the dwell 
time Td serves as the relevant long-time cut-off for quan- 
tum effects. (The TE-dependence of weak localization 
and conductance fluctuations in the presence of dephas- 
ing is considered in Refs. mii.) Our extension of the 
trajectory-based approach to arbitrary ballistic conduc- 
tors follows earlier work by Smilansky and coworkers 
who carried out a similar program for the diagonal ap- 
proximation to spectral fluctuations in closed quantum 
systems. Our final results are general expressions re- 
lating the Fano factor, weak localization, and conduc- 
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tance fluctuations to coarse-grained propagators of the 
classical dynamics in the ballistic conductor. We show 
that these general expressions reproduce known results 
for ballistic quantum dots and use our results to find the 
TE-dependence of the signatures of quantum transport in 
a quasi-one dimensional Lorentz gas. Our results for the 
Fano factor and the weak localization correction to the 
conductance agree with general expressions obtained in 
the field-theoretical formalism; The conductance fluctua- 
tions have not been calculated using the field-theoretical 
formalism, so that a comparison is not possible. 

Before we proceed with the exposition of the theory 
and a discussion of the results, two remarks about the 
specific form of the semiclassical limit and the appropri- 
ate classical propagators need to be made. First, we note 
that universality of quantum effects can be expected only 
in the limit Xp <C I, where I = vpT is the elastic mean 
free path, cf. Eq. ([T]) above. Since Af is proportional 
to Planck's constant h, this limit is equivalent to the 
semiclassical limit ?i — > 0. The Ehrenfest time, however, 
depends logarithmically on h through the ratio Ls/Xf, 
cf. Eq. (21) above. Since a nontrivial Ehrenfest-time de- 
pendence of quantum transport requires that te ~ td, 
one thus needs a semiclassical limit in which td grows 
logarithmically while sending ?i — > 0. An ?i-dependent 
increase of td requires an ?i-dependent change of the clas- 
sical dynamics. For the two examples considered here, a 
quantum dot and the Lorentz gas, this corresponds to an 
?i-dependent reduction of the size of openings or an in- 
crease of the system size L, respectively. Since the final 
results will be independent of the details of the classi- 
cal dynamics and since the required ?i-dependences are 
rather weak (only logarithmic in h), we believe that this 
slight modification of the classical dynamics while taking 
the semiclassical limit 7i ^ is inconsequential. 

The classical propagators appearing in our results will 
be coarse-grained both with respect to time and with 
respect to the phase space coordinates. The coarse- 
graining with respect to time is for a time window of 
order A~^; the coarse-graining with respect to the phase 
space coordinates corresponds to a distance Lg, which is 
the distance below which the classical dynamics can be 
linearized. The same distance Lg also appears in the def- 
inition of the Ehrenfest time. Since both and 
Lg/vp are much smaller than td, such coarse-grained 
classical propagators are sufficient for a description of 
quantum transport. The advantage of coarse-grained 
classical propagators is that no subtle phase space corre- 
lations need to be accounted for when evaluating the final 
expressions. In this respect, our final expressions differ 
from those in Refs. [13 and in which the weak local- 
ization correction and the Fano factor are expressed in 
terms of classical propagators in which quantum correla- 
tions are still implicit. (Vavilov and Larkin pointed out 
how these implicit correlations can be made explicit;'*^ 
the resulting theory has a form equivalent to the one 
presented here.) 

The remainder of this article is organized as follows. In 



Sec.|TT]we summarize the essential elements of the semi- 
classical formalism. In Sees. |IlT] and llVi we first discuss 
the Ehrenfest-time dependence of the Fano factor and 
the weak localization correction to the dc conductance. 
These sections show how the trajectory-based formalism 
is applied to ballistic conductors with an arbitrary geom- 
etry. In Sec.|V]we then turn to conductance fluctuations. 
The relation to the fluctuations of the density of states 
and the Gutzwiller trace formula is discussed in Sec. IVII 
We conclude in Sec. IVIII 

II. DEFINITION OF THE PROBLEM AND 
SEMICLASSICAL FORMALISM 

In our calculation, we consider a ballistic conductor 
coupled to electron reservoirs through two ballistic con- 
tacts. Transport is described using the scattering matrix 
S. Since there are two contacts, S has a block struc- 
ture S — Sjij, where the indices j' and j label the two 
contacts. 

Here and in the remainder of this article, primed variables 
refer to the exit contact. The dimension of the block Sj'j 
is Njr X Nj , where Nj is the number of channels in the 
jth contact, j = 1, 2; The dimension of the full scattering 
matrix S is N ^ Ni + N2. 

The conductance of the device is written 



where g is the dimensionless conductance, 

g = trS2iSl^. (5) 

We are interested in the interference correction (Sg) to 
the ensemble average {g), as well as the fluctuations of the 
conductance with respect to fluctuations of an external 
parameter that affects phases accumulated by the elec- 
trons, but not their classical dynamics. An example of 
such a parameter is the magnetic flux through an insulat- 
ing region in the sample interior or the Fermi energy for a 
device in which all scattering is from boundaries or scat- 
terers with a sharp potential profile. (The conductance 
fiuctuations with respect to variation of the classical dy- 
namics has been considered, e.g., in Rcfs. .l9.,43i .) We also 
consider the Fano factor for the shot noise power,— 

/ = —tr 5215*21 >S'22- (6) 

The Fano factor is self-averaging, and may be calculated 
by taking the ensemble averages of numerator and de- 
nominator separately. 

As discussed in the introduction, we calculate these 
observables in a semiclassical limit ?i — > in which the 
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FIG. 2: Hierarchy of time scales in two dimensional ballistic 
conductors in the semiclassical limit ?i — > 0. At finite tem- 
peratures or for time-dependent transport, the role of td is 
taken over by the minimum of td, the dephasing time r^, or 
the inverse frequency uj~^. 



ratio of the Ehrenfest time te and the dwell time td is 
kept fixed. Since rg grows logarithmically while sending 
h 0, this means that must also grow when the limit 
?i — *■ is taken. As a result, the relevant time scales 
can be grouped into four well-separated categories. The 
smallest time scale is Xf/vf, the microscopic quantum 
time scale of the problem. The next set of time scales 
consists of the classical time scales that do not grow log- 
arithmically with h in this semiclassical limit, such as 
the Lyapunov time or the elastic mean free time 
r. These separate the regime of the non-universal short- 
time electron dynamics and the universal long-time elec- 
tron dynamics that eventually determines the universal 
magnitude of the signatures of quantum transport. The 
third group of time scales consists of td and te. In the 
semiclassical limit taken here, these grow oc ln(l/?i) upon 
sending h ^ 0. The largest time scale is the Heisenberg 
time th = 27r?i/A, where A is the mean level spacing, 
which grows ex h~^. The results derived in the following 
three sections will be exact in the limit h ^ for this 
separation of time scales. It is the parametric separation 
between td and r that allows the use of coarse-grained 
classical propagators and removes the dependence on the 
non-universal short-time dynamics in the sample. 

One should note that, when taking the limit ?i at 
fixed te/td, the channel numbers A'^i and N2 and, hence, 
the conductance g diverge. (This divergence is linked to 
the hierarcy discussed in the previous paragraph since 
g ~ th/td-) The divergence of g ensures that effects 
related to Anderson localization can be ruled out. This 
divergence does not affect the quantities of interest to us, 
however, because the interference correction (Sg) to the 
average conductance, the variance of the conductance, 
and the Fano factor / remain finite. 

Our calculations are done using an expression of the 
scattering matrix Sj'j as a sum over classical trajecto- 
ries a that enter the sample through contacts j and exit 
through contact j' ^^^Sl 



j' j )m'm 



1/2 



Here m! and m label the propagating modes in the exit 
and entrance leads, respectively, and Wji and Wj are the 
widths of the entrance and exit contacts. The compo- 
nents p'^ and pj_,q of the momentum perpendicular to 



the lead axis are taken to be compatible with that of the 
modes m! and m in the corresponding leads. 



^ ±Trnm'/Wf, 
— ±TrUm/Wj, 



m! = 1, 
m — \, . 



(8) 



Further, Sa is the classical action of trajectory a and 
is its stability amplitude. The latter is defined as 



Aa — 



dy 



-1/2 



(9) 



where y is the coordinate perpendicular to the axis of 
the entrance contact and the partial derivative is taken 
at constant pj^. For simplicity of notation, the Maslov 
index and other phase shifts are included in the definition 

of Sa- 

Substituting the semiclassical expression for the scat- 
tering matrix ([7]) into Eqs. ([5]) and © one obtains semi- 
classical expressions for the conductance g and the Fano 
factor /. Since both Eqs. © and ^ contain products 
of the scattering matrix and its hermitian conjugate, the 
resulting expressions contain multi ple summations over 
classical trajectories. Following Ref.l45l we simplify these 
expressions in three steps. First, we note that trajectories 
that appear in neighboring factors of S and S*^ belong to 
the same transverse modes upon entry or exit, i.e., the 
magnitude of their transverse momenta is equal upon en- 
trance and/or exit. The case of opposite transverse mo- 
menta, however, is accompanied by a fast-varying phase 
factor which, if summed over, disappears in the semiclas- 
sical limit h Oi^i^ Hence, we only need to consider 
the case of equal transverse momenta. Second, in the 
limit h ~* 0, the summations over quantized transverse 
momenta p± and can be replaced by integrations. 
And third, locally the canonically conjugate coordinates 
p± and y can be replaced by the conjugate coordinates 
s and u, which are the stable and unstable phase space 
coordinates of the chaotic classical dynamics inside the 
conductor, defined for a Poincare surface of section at the 

lead opening. We choose their units to be equal, so that 

1/2 

both s and u have the same units as ?i ' . This coordinate 
transformation is accompanied by a Legendre transform 
of the classical action and the corresponding redefinition 
of the stability amplitude. The Legendre transformed 
action Sa is a function of the stable phase space coor- 
dinate Sa upon entrance and the unstable phase space 
coordinate u'^ upon exit. Its derivatives determine the 
remaining two coordinates upon entry and exit. 



dSa 

ds„ 



dSa 
du' 



(7) The stability amplitude reads 



A,. 



du' 



dUa 



-1/2 



(10) 



(11) 



where the derivative is taken at constant Sa- In this 
formulation of the semiclassical theory, trajectories a and 
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/3 belonging to factors Sij and 5*1; have equal stable phase 
space coordinates Sa = S/s in the entrance contact if j = 
I and Sij left-multiplies Slj, whereas they have equal 
unstable phase space coordinates — u'p in the exit 
contact if i = k and Sij right-multiplies Ski- We then 
arrive at the following semiclassical expression for g, 



N 

2TTh 



j ds^du'^^AaApe 

a,0 



(12) 



where the classical trajectories a and /3 run between con- 
tacts 1 and 2 with Sfs = Sq and u'^ — u'^. For the Fano 
factor / we find similarly 



/ = 



N 



dsadu'^dsjdu'^ 



X 2^ Ac,AfiA^Ase''-^"-^''+^-'-^'^/'', (13) 

a,/?, 7(5 

where the trajectories a and (3 connect contacts 1 and 2, 
the trajectories 7 and S connect contact 2 to itself, and 
the coordinates of the trajectories P and 6 satisfy Sfs = Sj, 
ss — Sa, u'p — u'^, and u'^ = u'^. These expressions will 
be the basis of the calculations of the next sections. 



III. SHOT NOISE 

In order to establish our methods and relevant approx- 
imations we first calculate the Fano factor /. 

A. Encounter in sample interior 

Technically, the simplest avenue to a semiclassical cal- 
culation of the Fano factor / is to use the unitarity of the 
scattering matrix to write / as 



/— -tr 52i>5'22'S'i2S'ii- 



(14) 



Since / is self averaging in the limit /i ^ 0, we may aver- 
age numerator and denominator in Eq. p4p separately. 
Using the semiclassical expression for S, we then find 



/ = - 



TV 



{27:hf{g) 



j dsadu'^dsjdu'^ 



AaApA^Ase'^^''-^''+^-''^''>/'', (15) 



where the trajectories a, (3, 7, and 5 connect the con- 
tacts 1 and 2, 2 and itself, 2 and 1, and 1 and itself, 
respectively. 

Only trajectories a, (3, 7, and 5 for which the total 
action difference A5 = Sa — + S^/ — Ss is of order 
fi systematically contribute to /. Such small action dif- 
ferences occur only if the trajectories a and 7, on the 
one hand, and /3 and (5, on the other hand, are piecewise 




FIG. 3: (Color online) Left: Schematic drawing of a set of four 
trajectories a, {3, 7, and 5 that contribute to the Fano factor 
/. Right: Detail of the encounter between the trajectories 
a and 7 together with definitions of the phase space points, 
coordinates, and partial encounter durations used in the main 
text. The true trajectories may involve specular refiection 
at the boundaries or at scatterers inside the sample. The 
separation between the trajectories is magnified for clarity; 
The true encounter regions involve trajectories separated by 
a sub-macroscopic distance which can not be resolved in a 
figure. 



identical, up to classical phase space distances of order 
Ti'^'^ or lessi ^°'^^ [Phase space distances are defined as 
max(|As|, |Au|), where s and u are the stable and un- 
stable phase space coordinates along a Poincare surface 
of section.] Since pairs of trajectories cannot be iden- 
tical for their entire duration because of the particular 
requirements on the entrance and exit contacts — see 
the text below Eq. ([TS]) — , this is possible only if the 
four trajectories undergo a 'small angle encounter', as 
shown schematically in Fig.[51 ^°d^'^^i'^^ In the encounter, 
the phase space distances between a pair of trajectories 
is of order ?i°c^~^° or less, where a > and c is a char- 
acteristic phase space distance below which the chaotic 
classical dynamics can be linearized. In the center of the 
encounter one has a < 1/2, while a = at the beginning 
and end of the encounter. All four classical trajectories 
are correlated for the entire encounter stretch because 
their phase space distance is below the cut-off c through- 
out the encounter. Before the encounter, the trajectories 
a and 5, and 7 and (3 are identical up to a quantum 
uncertainty, a phase space distance below ?i/c; After the 
encounter, a is paired with (3 and 7 is paired with (5, again 
with phase space distances up to a quantum uncertainty. 

Having identified the relevant sets of four classical tra- 
jectories contributing to (/), it remains to perform the 
trajectory sum in Eq. ([TS]). Since the trajectories (3 and 
5 are fully determined once a and 7 are specified, it is 
sufficient to sum over the relevant pairs of trajectories a 
and 7. The summation over such pairs of trajectories is 
performed by picking a reference point A on a somewhere 
along the encounter. In a two-dimensional ballistic con- 
ductor, one needs three classical coordinates to specify 
A (two position coordinates and the direction of propa- 
gation). Once A is chosen, the trajectory a is fixed. We 
fix 7 using the coordinates at which 7 passes through a 
Poincare surface of section taken at A. Following Refs. 

we parameterize a position on the Poincare 
surface of section using the stable and unstable coordi- 
nates of the classical dynamics at A. The phase space 
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point A is taken to be the origin of the coordinate sys- 
tem; The stable and unstable phase space coordinates of 
the point where 7 passes through the Poincare surface of 
section at A are labeled sa and ua, respectively. The pre- 
cise values of sa and ua will depend on where we choose 
the reference point A along the encounter. Moving A 
along the encounter, s and u change oc exp(±At), where 
A is the Lyapunov exponent of the classical dynamics in 
the sample. The encounter duration tcnc is defined as 
the amount of time during which the phase space dis- 
tance max(|s^|, I MA I) is less than the classical cut-off c, 
i.e., 

1 c 1 c 

tenc = ts+tu; ts^-ln- — -, t„ = - In -. (16) 

A \sa\ a \ua\ 

The action difference AS can also be expressed in terms 
of the stable and unstable phase space coordinates sa 
and UAr^ 



AS = saua- 



(17) 



Notice that neither tcnc nor AS depend on the choice of 
the reference point A. 

At this point, the summation over trajectories a and 7 
can be written as an integral over A, sa, and MAr^"^"^ 



1 



iAS/h 



(27r?i)2 to 



dA / dsAduA — p(A;sa,ua), 

(18) 

where the trajectory density p{A] sa, ua) is a product of 
delta functions that selects only those phase space points 
A and coordinates {sa, ua) for which the classical motion 
at A originates at contact 1 and ends at contact 2, while 
the classical motion at a point a phase space displacement 
(sA, Ua) away from A originates at contact 2 and ends at 
contact 1. The factor tone in the denominator cancels a 
spurious contribution to (/) from the freedom to choose 
A anywhere along the encounter. 

In a closed quantum system, phase space integrations 
similar to those of Eq. (fT8|) can be performed using the 
Hannay-Ozorio de Almeida sum rule.— Here, we take 
a different approach, following Refs. and lU, and re- 
place the exact trajectory density p by its 'statistical av- 
erage' {p), which is expressed in terms of classical prop- 
agation probabilities. The statistical average is taken 
with respect to small variations of the position in phase 
space^Siii and/or with respect to small fluctuations of 
the shape of the conductor. When replacing p by its sta- 
tistical average, we need to take into account that the 
propagation of the two trajectories a and 7 is correlated 
inside the encounter. In order to make these correla- 
tions explicit we introduce phase space points B and 
C on a at the beginning and ends of the encounter. 
The propagation times between the points B and A 
and A and C are ts and tu, respectively. The trajec- 
tory 7 passes through Poincare surfaces of section at B 
and C at coordinates {sb,ub) = [sa^^*' ,UAe-~^*') and 
{sc,uc) = (sAe~'^*",UAe'*'*"), respectively. Beforehand 



after C the trajectories a and 7 are uncorrelated, so that 
the probabilities that a and 7 originate from/end at the 
appropriate contacts factorize. Hence we have 

{p{A;sa,ua)) = j dBdCP{A,B-ts)P{C,A-tu) (19) 
X P{B, 1)P{B*,2)P{2, C)P{1, C*). 

Here P{A, B;ts) is the probability that a classical tra- 
jectory starting in the phase space point B reaches A 
after a time ts, P{C, A; tu) is the probability that a clas- 
sical trajectory starting in A reaches C after a time tu, 
P{B,j) is the probability that a trajectory at B origi- 
nates from contact j, and P{j', C) is the probability that 
a trajectory at phase space point C ends at contact j' . 
(Here and in the remainder of this article we use a semi- 
colon ";" to separate phase space and time arguments 
of the classical probabilities. The absence of a time ar- 
gument indicates that the classical probability has been 
integrated over time.) The phase space points B* and 
C* are a phase space displacement {sb,ub) and (sp, uc) 
away from B and C, respectively. For the semiclassical 
limit taken here (dwell time td larger than the Lyapunov 
time A^^ by a factor logarithmically large as 7i ^ 0), a 
phase space displacement over a distance ~ c does not 
affect propagation probabilities, so that we can coarse- 
grain the probabilities P{B,j) and P{j' , C) over a phase 
space volume of size ~ c^/A and neglect the difference 
between B and B* or C and C*. 

We eliminate the phase space coordinates sa and ua 
in favor of phase space coordinates sb, ub at a Poincare 
surface of section taken at B. This is done with the help 
of the variable chang o^^'^^ 



sa — ever, ua = cx/va, 



(20) 



where a — ±1, \x\ < 1, and |a;| < d < 1. In terms 
of the new integration variables one then has sb — ca, 
Ub — cx/a. Further, the action difference AS — (?x and 
the encounter duration 



tone = A-iln(l/|x|). 



(21) 



Now, the initial reference point A can be integrated out 
using 

y dAP{A, B; ts)PiC, A; tu) - P{C, B; W). (22) 

The integral over v can be done as well and cancels the 
factor tcnc in the denominator of Eq. ^TE\\ . We then find 



{2nh)^g) 



J dBdC J dxY^ e"''' 



/h 



X P(C, B; te„c)P(S, 1)P{B, 2)P{2, C)P{1, C). 

Next, wc sum over a and perform a partial integration 
to x, with the result 

(/) = -^^2^ / dBdCP{B,l)P{B,2)P{2,C)P{l,C) 

(23) 



^ dx . xc^ d 



JO X n CEonc 
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The X integration in Eq. ([23)1 converges for x ^ Ti/c^. 
The probabihty P{C,B; <enc) depends on x through the 
ratio icnc/TD only, where td is the typical dwell time. We 
write icnc = te + where 



te = A~iln(cV?i) 



(24) 



is the Ehrenfest time and the remainder St is of order A~^. 
Since we take the semiclassical limit ?i — s- at fixed ratio 
te/td, one has St/rD — > 0, so that <onc may be replaced 
by Te in the argument of P{C, B;tcnc)- We then arrive 
at the final result 



^JdBdCPiB,l)PiB,2) 



X P(2, C)P{1, C)—P{C, B- te). (25) 
Cte 



The factor dP{C , B\t-e) / Bt-e is the ballistic counterpart 
of the 'Hikami box' in the diagrammatic perturbation 
theory of disordered conductors fiSii^ 

This expression for the Fano factor is equivalent to 
a similar expression obtained in Ref. 15 using the field- 
theoretic approach. This need not be obvious at first 
sight, because the final expression in Ref. 15 involves an 
integration over a single phase space point A in the sam- 
ple interior only. In the language employed here, this 
phase space point is located at the center of the en- 
counter, where the distance between the trajectories a 

1/2 

and 7 is of order h ' . The classical propagators in Ref. 
[Tsl are not coarse grained, but they should be evaluated 
in the presence of the diffraction off an additional weak 
potential, which effectively amounts to a smearing of the 

1/2 

classical propagators over a phase space distance ^ h ' , 
see Ref. |42. Taking into account that classical trajecto- 
ries are still strongly correlated at a phase space distance 
~ h^^^ and that it takes a propagation time te/2 away 
from A for such correlations to disappear, one arrives at 
a structure similar to our Eq. (|25p above. 



B. Encounter touching the lead opening 

Alternatively, the Fano factor could have been calcu- 
lated from the original expression ([6]) or its semiclassical 
version, Eq. p^ . That the semiclassical evaluation of / 
from a different expression than the one used in the cal- 
culation above gives the same result was shown in Refs. 
[23II47I and m for the case of a ballistic quantum dot. We 
now extend this verification to the general case. 

If / is calculated from Eq. (fT5|) . the trajectories a and /3 
enter the sample through contact 1 and exit through con- 
tact 2, whereas 7 and S enter and exit through contact 2. 
The fact that all trajectories exit the sample through the 
same contact leads to an additional contribution that did 
not exist in the calculation described above. This addi- 
tional contribution arises from small-angle encounters of 
the four trajectories a, /3, 7, and S where the encounter 
touches the exit contact, i.e., the phase-space distance 



between the trajectories a and 7 is less than c upon exit 
from the sample. (The reason why such a contribution 
does not exist for the calculation described above is that 
the trajectories a and 7 exit through different contacts 
in that case.) 
Hence, we write 



(26) 



where (/)^^^ denotes the contribution from encounters 
that do not touch the contacts and (/)^^^ denotes the 
contribution from encounters that touch the exit contact. 
Proceeding as before, we find 



(1) 



J dBdCP{B,l)P{B,2) 



xP{2,Cf—P{C,B;TE). 



(27) 



In order to calculate (/) one again proceeds by select- 
ing a phase space point A at an arbitrary point during 
the encounter, as well as a Poincare surface of section 
with stable and unstable phase space coordinates sa and 
ua- The propagation time i„ between A and the lead 
opening is bounded by 



1 , c 
tu< ^In -. r 

A \ua\ 



(28) 



whereas the propagation time between A and the begin- 
ning of the encounter is tg = A^^ln|c/s^|. The total 



action difference AS = Sa 
before.''^ We thus find 



Sp + S~f — Ss 



saua, as 



If /"C pis 

if)^^^ ^ j^JdAjjsAduAj^ 



-p{A] sa,ua), 
(29) 



where <onc — ts+tu and p{A; s^, ua) is a product of delta 
functions that selects only the appropriate phase space 
points A and coordinates sa and ua- We again replace 
p by its statistical average. Defining a point B at the 
beginning of the encounter (the end of the encounter is 
at contact 2), we have 

{p{A-s,u)) = j dB j dtuP{2,A;t^) 

X P{A,B:U)P{B,1)P{B,2), (30) 

where the integration domain for tu is given by Eq. (|28|) 
above. We change variables 

UA = cx/va, SA^cva, t„ = tone - A"Mn(l/i;), (31) 

with a = ±1, < 1, and \x\ < v < 1, so that the 
phase space coordinates for a Poincare surface of section 
at point B are ub = cx/a and = crc. Then we inte- 
grate over the phase space point A and over the phase 
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space coordinate v and sum over a 
2c^X 



if) 



(2) 



{2nhr{g) 

r-A"^ ln\l/x 



J dBP{B,l)P{B,2) J dxe""'/'' 



(^^cnc-f (2, -S, ^cnc)- 



(32) 



Performing a partial integration to x, and then integrat- 
ing to X, we find the resuh 

^^^^'^ ^ / 2)^(2' 

The sum (/> = (/>(i) + agrees with the resuh ^ 

obtained previously, as one easily verifies using current 
conservation, 



P(1,C) + P(2,C) = 1, 



(34) 



and the equality 



J dCdrPU, C)P{C, B- t) + P(j, B- t) =0, j = 1, 2. 

(35) 



C. Quantum dot 

In a quantum dot, one has 

P(j,C) =iV,/7V, J = 1,2, 

and 

P(A,P;r) = le-^/^-, 



(36) 



(37) 



where f2 = 2TThNTD is the dot's total phase space vol- 
ume. Also, to leading order in A^i and N2, the average 
conductance is the series conductance of the two point 
contacts, i.e., 



{9) = N,N2/N, 
independent of the ratio te/td. Hence 



(38) 



(39) 



in agreement with the Fano factor previously obtained in 
Refs. [H andlH. In the limit te/td 0, Eq. jSH) sim- 
plifies to the Fano factor obtained from random matrix 
theoryi^ 



D. Lorentz gas 

We now apply the general expression (|25|) to a quasi- 
one dimensional Lorentz g random collection of disc- 
like scatterers placed in a ballistic electron gas with a 
width W much smaller than its length L. In a Lorentz 




FIG. 4: Ehrenfest-time dependence of the Fano factor / in 
quasi-one-dimensional Lorentz gas. 



gas, the coarse-grained electron motion is diffusive and 
described by a diffusion constant D. The typical dwell 
time for electrons transmitted through the Lorentz gas is 
taken to be 



(40) 



[The factor tt^ in Eq. (l40l) is added for notational conve- 
nience.] Taking the limit ?i at a constant ratio te/td 
and constant sample width W implies that L scales pro- 
portional to [ln(l/?i)]^/^ if ?i ^ 0. The relaxation of the 
momentum and the transverse coordinate y of an elec- 
tron occur on time scales ~ and ~ W'^ /D, respec- 
tively, which do not scale with h. Since the relaxation 
of X is parametrically slower, we may coarse- grain the 
propagation probabilities appearing in Eq. (|25p over the 
propagation angle <j) and the transverse coordinate y. 

In order to evaluate Eq. ([^5]) , we write the phase-space 
elements dB and dC as 

dB — 2nh^i'dxBdyBd(j)B, dC = 2TTh^vdxcdycd(j)c, 

(41) 

where the angle (j) denotes the propagation angle and u is 
the (two-dimensional) density of states. Integrating over 
(t'B, yCi and (j>c, one then finds 



2TThiyW 



j dxdx'Pix,l)P{x,2) 



X j^P{x', X- te)P(1, a;')P(2, x'), (42) 

OTE 



where W is the wire width, P{x,x';t) is the diffusion 
propagator in one dimension, and P(j, x) is the proba- 
bility that diffusive motion originating at position x exits 
the sample through contact j, j = 1.2. For a sample of 
length L one has 



P{x,x';t)^-e{t)J2 



-u^t/rn ■ A^""^ • Mi's: 

e ^ ^'"^^ sm - — sm 



(43) 
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where td was defined in Eq. (|40p above and 9{t) = 1 
a t > and otherwise. The probabilities to escape 
through the contacts 1 and 2 read 



P(x,l) = P{l,x) = 
P(x,2) = P{2,x)^ 



L — X 



so that 



Since 



(/) = 



2T:hvDW 



32 



L{9) 



(g) 



fi odd 



2TrhiyDW 



L 



(44) 
(45) 

(46) 
(47) 



for a two-dimensional wire of length L and width W , we 
find that the Fano factor reads 



^ 32 

fi odd 



(48) 



The TE-dependence of the Fano factor is shown in Fig. 
m In the hmit te/td 0, Eq. agrees with the well- 
known result _F = 1/3 for a disordered wirci^ For finite 
te, the Fano factor is smaller, but the dependence on te 
is not described by a single exponent, as in the case of a 
chaotic quantum dot. 



IV. WEAK LOCALIZATION 

The semiclassical calculation of the weak localization 
correction (Sg) to the ensemble averaged conductance (g) 
is essentially equal to that of the Fano factor for the shot 
noise power. Starting point is the semiclassical expres- 
sion (|12p . The relevant trajectories contributing to (Sg) 
in the semiclassical limit are shown in Fig. [5l They are 
pairs of trajectories a and (3 with a small-angle intersec- 
tion and a loop, such that a and /? are identical (up to 
quantum uncertainties) before and after the intersection, 
but traverse the loop in opposite direct ions J^i^S 

Although the summation is over pairs of classical tra- 
jectories, the trajectory /3 is fully specified once a is fixed. 
Hence, the calculation of (Sg) requires a single summa- 
tion over all relevant trajectories a with small-angle self 
intersections. For chaotic cavities, such sums were first 
performed by Sieber and Richteri^^i^ Here we extend 
this calculation to arbitrary geometries. As in the previ- 
ous section, the only assumption of our calculation is that 
the motion is locally chaotic with Lyapunov exponent A. 

Proceeding as in the previous section, we parameter- 
ize a using a phase space point A taken at the first 
passage along the small-angle encounter. We take a 
Poincare surface of section at A, and label the unstable 
and stable phase space coordinates of the second pas- 
sage of the time-reversed of a through the encounter as 




FIG. 5: Left: Schematic drawing of a pair of trajectories a 
and /3S that contributes to the weak localization correction to 
the conductance. Right: Detail of the self-encounter of one 
of a together with the definition of the phase space points A, 
B, and C used in the main text. 



{satUa)- The encounter duration fonc aud action differ- 
ence A5 = Sa — Sp are given by Eqs. (fTB|) and (fTT)) of 
the previous section. With this, the summation over tra- 
jectories a can be written as an integral over A^ sa, and 



UAf 



59 



dA / dsAduA 



-p{A;sa,ua), (49) 



where the trajectory density p{A;sa,ua) selects only 
those phase space points A and coordinates sa, ua for 
which the classical motion from A originates from the 
contact 1, continues to a point a phase space displace- 
ment {sa,ua) away from the time-reversed of A, and 
ends at contact 2. The factor ienc in the denominator can- 
cels a spurious contribution from the freedom to choose 
A anywhere along the encounter. 

As in the previous section, we replace the exact tra- 
jectory density p by its statistical average (p). We again 
introduce the encounter beginning B and its end C. Clas- 
sical propagation probabilities outside the encounter fac- 
torize, so that 

{piA;s,u)) = P{A,B,ts)PiC,A;tu) 

X P{2,B)P{B,1) J dtP(C,C;t), (50) 

where B and C denote the time-reversed of B and C, 
respectively, up to phase space displacements of order of 
the cut-off c. As before, we used probability distributions 
that were coarse-grained over a phase space volume ~ 
c^/ A. Repeating the variable changes of Eq. ((20|) and the 
manipulations following that equation, we arrive at the 
result 

{6g) = J dBdCP{2,B)P{B,l) 

X J dtP(C,C;t)-^P{C,B;TE). (51) 

This expression agrees with the result previously de- 
rived by Aleiner and Larkin for an arbitrary ballistic 
conductor fiS, see also Refs. [23. 
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The factor P{C, C; t) is the equivalent of the 'Cooperon 
propagator' in the theory of weak locahzation in disor- 
dered conductors. In the presence of a magnetic field, 
P{C,C;t) should be multiplied with the phase factor 
exp(2ie$/ft,c), where $ is the flux enclosed in the loop 
connecting C and C. If $ ^ hc/e for a typical loop, the 
weak localization correction is suppressed. 

Instead of calculating the interference correction (Sg) 
from the total transmission, as was done above, one may 
also calculate {6g) from the reflection coefflcient, 



one flnds 



= iVi -trS'iiS'li. 



(52) 



If one proceeds this way, there are two interference con- 
tributions to (Sg), 



(Sg) ^ {Sgy'^ + {Sgy'K 



(53) 



These two quantum corrections are known as the 
weak-localization correction to reflection and coherent 
backscatteringi^Si^i They give a positive and negative 
quantum correction to g, respectively. The flrst contribu- 
tion arises from trajectories that have a small-angle self 
encounter that fully resides in the sample. It has a form 
similar to Eq. (|5ip above, 



(Sg) 



(1) 



dBdCP{l,B)P{B,l) 
dtP{C,C-t)^P{C,B-T^). (54) 



The second contribution is from trajectories that have a 
small-angle self encounter that touches the contact. Its 
calculation proceeds analogous to the derivation of Eq. 
and has the result 



(2) 



dC dtPiC,C;t)P{C,l;TE). (55) 



(with C replaced by B) one 
(^g)^^) agrees with the result 



Using Eqs. ^ and ^ 
verifles that (Sg) = (^g)*-^^- 
([51]) obtained previously i^°i^^ 
Applying Eq. ([STjl to a chaotic quantum dot, one sets 



P{C, B- 1) = P{C, C;t) = ^ exp(-i/rD), (56) 

where f2 — 27r?iA^TD is the total phase space volume of the 
quantum dot, and P(j, B) = P{B,j) = Nj/N, j = 1,2. 
The factor fl cancels from the expression for {6g) , so that 



(-5.9) 



N1N2 



-te/td 



(57) 



in agreement with Refs. 22, gHH (see also Refs. . In 

the limit te ^ td Eq. (j57p reduces to the weak localiza- 
tion correction obtained within random matrix theory.'^ 
The weak localization correction (Sg) for a quasi-one 
dimensional Lorentz gas of length L and classical diffu- 
sion constant D is calculated by by replacing the coarse- 
grained probabilities appearing in Eq. (|5ip by propaga- 
tors of a one-dimensional diffusion process. In that case, 



{6g) = - 



E 

fj.. odd 



32 



(58) 



In the limit te ^ / D of short Ehrenfest times, this 
result simplifies to the well-known weak localization cor- 
rection SG — — (l/3)(2e^//i) of diffusive quantum wiresj^ 
If Te and /D are comparable, however, weak localiza- 
tion is suppressed. The suppression has the same func- 
tional form as the suppression of the Fano factor for the 
shot noise power, see Fig. |4l 



V. CONDUCTANCE FLUCTUATIONS 

The same framework can be used to calculate the con- 
ductance fiuctuations in arbitrary ballistic conductors 
in the semiclassical limit. Conductance fluctuations are 
characterized by the correlation function 



Kit) 
K{e~e') 



1 

2^ 



dLuK{hLu)e" 



{g{e)g{e')) - {g{s))' 



(59) 
(60) 



The correlation function K{t) determines the variance of 
the conductance at flnite temperatures. 



varg(r) = / dt 



Kit). 



(61) 



At a flnite temperature, dephasing from electron-electron 
interactions further suppresses the conductance fluctua- 
tions. The effect of thermal smearing considered here is 
dominant, however, since typically Tt^ 3> fi.^ (See Ref. 
28, for a discussion of the TE-dependence of conductance 
fluctuations in the presence of dephasing.) The Fourier 
transformed correlation function Kiuj) describes the en- 
ergy dependence of the conductance fiuctuations at zero 
temperature. 



A. Encounters in sample interior 

Following Refs. |2^ and[27l. where Kit) was calculated 
for a ballistic quantum dot at finite Ehrenfest time, we 
calculate the conductance autocorrelation function Kit) 
by setting 



gie) = Ni ~ tY Snie)Snie)l 
gie') - N2~tYS22ie')S22ie')^ 



(62) 



in Eq. (|5D|) . This formulation avoids the necessity of deal- 
ing with encounters that touch the contacts. (The case of 
encounters that touch the lead openings will be discussed 
later.) 

Using semiclassical expressions for 5*11 and S22, Kit) 
is expressed in terms of a quadruple sum over classical 
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c 



FIG. 6: (Color online) Schematic drawing of sets of four tra- 
jectories that contribute to the conductance fluctuations. In 
panels b and c the periodic reference orbit 7 is shown thick. 
The period t-, of y is t. 

trajectories and f3i, i = 1,2, each connecting contact i 
to itself. The contribution of each set of four trajectories 
involves the action difference 

AS = ~ 5/3^ + - 5^2- (63) 

Only combinations of four trajectories for which the to- 
tal action difference AS is of order h systematically con- 
tribute to vai g. 

The action difference AS is a function of the energies 
e and s' , 

AS{e,e') = AS{O,O)+e{Tc,,~Ti3,)+e'iTa,~T0,), (64) 

where Tq. and Tp. are the durations of the correspond- 
ing trajectories, i = 1,2. As we'll show below, for all 
trajectories that contribute to K{t) one has 

SO that AS is a function of the difference e ~ e' only. 
Moreover, because of the Fourier transform to e — e', cf. 
Eq. ([55)1 above, only trajectories with 

t = Tp^- Tqj = Ta^ - Tfj^ (66) 

contribute to K{t). 



There are three classes of trajectories that meet these 
criteria and contribute to K{t) in the semiclassical limit 
?i — > 0.— They are shown schematically in Fig. [6] for 
the case of positive t. (Trajectories for negative t can 
be found by interchanging aj and j3j, j — 1,2.) Figure 
shows four trajectories that undergo two successive 
encounters, separated by stretches of independent prop- 
agation of duration r and t + t, respectively. Figures 
and c show configurations of trajectories in which the tra- 
jectories ai and f3i differ by a periodic reference orbit 7 
of period r-y = t, which is in (3i but not in ai. The same 
periodic orbit is also the difference of (32 and 0:2 ■ The 
configuration of Fig.[6j;, which involves a single encounter 
of all four trajectories that extends before and after the 
encounter with 7, can be considered as originating from 
that of Fig. [5^ if the two encounters there were to overlap, 
i.e., if r < 0. No correlations exist away from the periodic 
reference orbit 7 for the trajectory configuration of Fig. 
|BJd. The two encounters in Figs. [5^ and b have a duration 
te; The single encounter in Fig.[BJ; has a variable duration 
ranging from te to 2rE. The trajectories of Figs. [5^ and 
b have their counterpart in the diagrammatic theory of 
conductance fluctuations in disordered metals j^^i^^'^'^i^^ 
the configuration of Fig. [Bj; exists in ballistic conductors 
at finite te only. In the presence of time-reversal sym- 
metry, three additional contributions to the conductance 
autocorrelation function appear, which are obtained by 
time-reversing the trajectories 02 and P2 in Figs. [6^-c. 
Configurations of interfering trajectories with more than 
two encounters give contributions to var G smaller by a 
power of A'^ and need not be considered in the limit h ^ 
at fixed t^/tj:) we consider herci^ 

In order to establish that the relevant configurations 
of trajectories indeed are the three configurations shown 
in Fig. [5] we divide the possible trajectory configurations 
into two groups: those that involve one or more revolu- 
tions around a periodic orbit (as in Figs, and c) and 
those that do not (as in Fig. [B^). Although this proce- 
dure ignores the close relation between the configurations 
of Figs. [B^ and c, it proves to be better suited to an unbi- 
ased evaluation of all possible trajectory configurations. — 
Calculation of the contribution of interfering trajectories 
without revolutions around a periodic orbit is straight- 
forward since the quantum interference correction from 
non-overlapping encounters factorizes.^^ The case of a 
single encounter was considered in Sec. IIIII Thus, we 
find that the contribution K{tY°''> of trajectories of the 
type shown in Fig. [6^ reads 



Kit)'-''^ = J dAdB j dCdDP{A,l)P{A,2) 
X j dTP[C,B-T + t)P{C,B;T) 



—P{B,A;te) 
^P{D,C;te) 



P{l,D)P{2,D), 



(67) 
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FIG. 7: (Color online) Schematic drawing of two of the interfering trajectories (one of each pair) for all possible contributions 
to K{t). The figure also contains the definitions of the phase space points (represented by solid dots), Poincare surfaces of 
section (thick arrows) and time intervals used in the main text (thin arrows). 

where we use the convention that P{-, •;t) = if r < 0. The positions of the phase space points A, B, C, and D, as 
well as of the four Poincare surfaces of section used in a calculation of Eq. (pT]) , are indicated in Fig. [7^. 

In order to describe the contribution from trajectory configurations that involve a revolution around a periodic 
orbit we use the trajectories ai and (32 and the periodic orbit 7 as a reference, see Figs. [6)d and c. The remaining 
trajectories f3i and a2 are fixed by the boundary conditions that these must be paired with ai and (32 before and 
after the encounter with 7, respectively, and that Pi and a2 have one extra revolution around 7. The calculation of 
the contribution from trajectory configurations that involve revolutions around a periodic orbit is technically more 
involved than that of the contribution of Fig. [6^ because of the possibility that the encounter of the trajectories ai 
and (32 with the periodic orbit 7 winds around 7 and, hence, overlaps with itself. In order to account for the possibility 
of overlapping encounters, we start from phase space points Ai and A2 taken at the periodic reference orbit 7 at the 
beginning of the encounters of the periodic orbit and the trajectories ai and (32 with the periodic center trajectory, 
see Fig. [T] At each of these phase space points we draw a Poincare surface of section, for which we use phase space 
coordinates (u;, Si), i = 1,2, referring to the unstable and stable phase space coordinates at which ai and (^2 pierce 
these Poincare surfaces of section, respectively. (The origins of the coordinate systems are taken at Ai and A2.) Since 
Ai and A2 are at the beginning of the encounter, we have Si = coi, with ai = ±1, i = 1, 2. We write Ui = cxij Oi. The 
trajectories a2 and (3\ pass through these Poincare surfaces of section at least twice and have phase space coordinates 
{uie~^*,Si) and (ui,Sie^^*), i — 1,2, upon their first and second encounter, respectively. The action difference then 
reads2^ 

AS = (w2S2-MiSi)(l-e-^*)-t(£-e'). (68) 

We are interested in times t ^ A^^ for which we may neglect exp(— Ai) in comparison to unity. 

We count periodic trajectories that consist of several revolutions of one shorter trajectory as separate trajectories. 
This correctly takes into account the contribution from interfering trajectories where the difference between interfering 
trajectories is more than one revolution around a periodic trajectory. 

Now we are ready to calculate the contribution of all trajectories that involve one or more revolutions around a 
periodic orbit. Repeating the steps used for the calculation of the Fano factor and the weak localization correction, 
we find 

K{tY'+'^ - (f^) dxidx2e^='(^i-^^)/'' ^ f dA,dA2p{Ai;A2;<7,,Xi;a2,X2;\t\). (69) 
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The trajectory density p selects only those trajectories for which the phase space points Ai and A2 lie on the same 
periodic trajectory 7 with period and for which the trajectories ai and [32 (which are determined by Ai, cti, and 
xi, and A2, (72, and X2^ respectively), enter and exit through contacts 1 and 2, respectively. 

As before, we replace the exact trajectory density p by its statistical average (p) and express (p) in terms of classical 
propagation probabilities. Hereto, we need the time T12 of propagation from A2 to Ai. We require |ri2| < |t|/2, thus 
allowing for negative T12 if the propagation time from Ai to A2 is shorter than the propagation time from A2 to Ai. 
We also introduce points Bi and B2 on the periodic trajectory, located at the end of the encounters of 7 with ai and 
/32, see Fig. [71 The travel time between Ai and Bi is 

isA,. = A-Mn(l/|x,|). (70) 

The order in which the four phase space points Ai, A2, -Bi, and B2 appear on the periodic trajectory 7 is determined 
by the times \t\, tu, tAB,i, and tBA,2- 

If the trajectories ai and (32 are already correlated at the point where they first meet the periodic trajectory, i.e., 
if the encounter of ai and (32 begins before the encounter of these trajectories with the periodic trajectory, we also 
introduce the phase space point C at the beginning of their encounter. This case is shown in Figs. and d. In that 
case the phase space points Ai and A2 are classically close, because the phase-space separation between ai and (^2 is 
sub-macroscopic if their encounter has begun before the encounter with 7. The travel time between C and Ai is 

r^c = max(0,A-Mn|l/s|), (71) 

where cs is the difference between the stable phase space coordinates of ai and P2 taken at the point Ai, 

s = CTi - cr2e~^^^^ (72) 

Similarly, if the encounter of ai and (32 extends beyond the encounter of these trajectories with the periodic trajectory, 
we also introduce the phase space point D at the end of that encounter. In that case, the phase space points Bi and 
B2 are classically close. The travel time between Bi and D is denoted tfb- 

When expressing (p) in terms of classical propagation probabilities, we need to distinguish between the cases with 
and without correlations between the trajectories ai and P2 while away from the periodic trajectory. Hereto we write 

(P> = {P)^'^ + {P)^'^ + {P)^'^ + {P)^'^- (73) 

The four terms correspond to the four cases shown schematically in Fig. [7jD-e. The first term (p)*^''' describes the case 
in which the trajectories ai and (^2 are not correlated away from the periodic trajectory 7. This is the case shown in 
Figs. [Tjs or[6}3. One finds 



(p)('') = dBidB2 dTi2P(Ai,l)P(A2,2)P(l,Si)P(2,B2)P^(Ai,A2,Bi,S2;|t|,Ti2,iBA,l,isA,2), 
J "'-|t|/2 

(74) 



where denotes the probability that the four phase space points Ai, A2, Pi, and P2 are on a single periodic 
trajectory 7 with travel times between these points as specified by the times \t\, tyi, ^ba,!, and Iba.i- (Note that 
^SA.i and tBA.2 may be larger than the period t of 7. In this case, the encounters of ol\ and (32 with 7 are 'wrapped 
around' the periodic reference trajectory 7.) The third term (p)^''-' is a correction for the case that the two trajectories 
are correlated before they reach the periodic trajectory. This case is shown in Fig. [Tjl. The average trajectory density 
(p)^'') reads 

(p)^"") = J dB,dB25(Ai - A2) I dn2 I dCP{C,l)P{C,2)P{AuC;tAc)~P{Ai,l)P{Ai,2) 

X P(l, Pi)P(2, B2)P-y{Ai,Bi,B2; \t\,tBA,i, tBA.2))- (75) 

The integrand implicitly depends on T12 through the time tac of propagation between C and Ai, cf. Eqs. (j7ip and 
([7^ . The integration domain for T12 is limited to those T12 for which tac appreciably differs from zero. This occurs 
if (Tif72 = 1 only, see Eqs. (|7T|) and ([7^ . In that case the integration range for T12 is cut off at |ti2| ^ . Since the 
Lyapunov time A""'^ is parametrically smaller than |t|, tBA,i, and tBA,2, we could identify Ai and A2 in the r.h.s. of 
Eq. (j75p and remove T12 as an argument of P-y. Similarly, one has 

(p)(^) = J dB J dTuP{Ai,l)P{A2,2) j dDP{D,B)P{l,D)P{2,D)~ P{\,B)P{2,B) 

X P^{Ai,A2,B-\t\,tBAA,tBA,2) (76) 
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and 



J dB6{Ai~A2) J dri2 J dCP{C,l)P{C\2)P{A,,C;tAc) - P{Ai,l)P{A,,2) 
J dDP{D, B)P{1, D)P{2, D) - P(l, B)P(2, B) 



P-,{Ai,B-\tltBA). 



(77) 



see Figs. [7h and e. The corresponding four contributions to K{t) need to be considered separately. 

For the first contribution Kit^^^ we first perform a partial integration to xi and to X2- The remaining integral over 
xi and X2 is convergent and sets tAB.i and tAB,2 equal to the Ehrenfest time te, up to corrections of order that 
are not relevant here. Hence, we obtain the result 



ib) 



J dAidA^dBM 'dTi2P(Ai,l)P(A2,2)P(l,Bi)P(2,B2) 



d 



d 



dtABA dtAB,2 



Py{Ai,A2,Bi,B2; \t\,Ti2,tAB,l,tAB,2) 



(78) 



tAB,l—tAB,2—TE 



where we chose the integration domain for ri2 to be < T12 < \t\ instead of the original choice — |t|/2 < T12 < \t\/2. 
For second contribution from (p)'"'-' one again first performs partial integrations to xi and X2. As in Eq. (|78p . the 
remaining integration over T12 is regular. However, the integration domain for the T12 integration is of order A~^, 
not \t\, so that the magnitude of this contribution to K{t) is a factor ^ l/'^l^l ^ 1 smaller than K{t)'^^\ Similar 
arguments apply to the contribution from {pY'^^ ■ The contribution K{t)^'^'> from (p)'-'^-' is not small, however, because 
the T12 integration is singular in this case. For the encounter of ai and (32 to extend before and after the encounter 
with 7, tAB.i and tAB,2 must be equal, up to multiples of \t\. Only the case Iab.i — tAB.2 needs to be considered, 
however, because only in this case the T12 integration is singular 1^ Also, one must have oi = (T2 and signui = signu2 
for the trajectories ai and f32 to be correlated both before and after their encounter with 7. We shift to integration 
variables 



X = wiai(l-e-^"i=)/c, tBA = A-iln(c/|Mi|), y ^ {ui ~ U2e^^^^) / c, 



(79) 



so that 



- ^[^)' I dxdy I dtBAC-'^^^y^/' 



J dAdB J dCP{C, 1)P(C, 2)P{A, C; Iac) - P{A, l)P{A, 2) 
j dDP{D, B; tDB)P{l, D)P{2, D) - P(l, P)P(2, B) 



P^{A,B-,\tltBA), 



(80) 



where we included two factors 2 for the two possible values of ui — 02 and the two possible values of signui. In the 
calculation of the action difference, we neglected x and y with respect to unity, which is allowed because the final 
integration converges for x, y ^ Ti/c^ <^ 1. The times Iac and t^B read 



tAC = X-^\n{ll\x\)-tBA, tDB=X'^H^I\y\)-tBA. 



(81) 



Again, we perform partial integrations to x and y. After partial integration, the integrals over x and y converge and 
we find tsc = ^Ac + tBA = te and tjjA = toB + tBA = te, up to corrections of order A""'^ that can be neglected when 
appearing as a time-argument in the classical probabilities. Hence 



K{t) 



J dtBA j dAdBdCdDP{C,l)P{C,2)P{l, D)P{2, D) 
X T^P(A, C; TE - tBA)^P{D, B; te - tBA)P^{A, B; \t\,tBA)- 

OTE OTE 



It remains to write down explicit expressions for the two probabilities P^ appearing in Eqs. ([75]) and 
probability Pj{A,B; \t\,tBA) that appears in Eq. ([5^ is given by 



(82) 
The 



P^{AB;\tltBA) ^ PiB,A;T)P{A,B;\t\-T), 



(83) 
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with 

T^tAB mod (84) 
As before, we use the conventfon that P(-, •;t) = if r < 0. Similarly, for the probability P-y of Eq. (|78|) we find 

P^(Ai, A2, Pi, P2; ri2, TAB, 1, TAB, 2) 
= P(Pi, ^1; ri)P(A2, Pi; T12 - Ti)P(P2, A2; T2)P(Ai, P2; \t\ - T12 - T2) 

+ P(yl2, ^1; Ti2)P(Pi, v42; Ti - Ti2)P(P2, Pi ; r2 + T12 - ri)P(Ai, P2; \t\ - T12 - T2) 

+ P(P2,^i,Ti2 +T2 - |t|)P(Pi,P2;Ti + |i | - ri2 - )P(yl2 , Pi , T12 - Ti )P(yli , ^2 ; |i | " ri2 ) 

+ P(A2, ^1, ri2)P(P2, A2; r2)P(Pi, P2; Ti - r2 - ti2)P(Ai, Pi; \t\ - n) 

+ P(Pi, ^1; Ti)P(P2, Pl, ri2 +T2- \t\ - Ti)P(A2, ^2; 1^1 " T2)P(Ai, ^2; \t\ " T^) 

+ P(P2,^l,Ti2+T2- |t|)P(A2,P2;|i| -T2)P(Pi,A2;Ti -ri2)P(Al,Pi;|t|-Ti) , (85) 

where the six terms correspond to the six permutations of the four phase space points Ai, A2, Pi, and P2 along the 
periodic reference trajectory 7 and 

Tj = tAB,j mod |t|. (86) 
The full conductance autocorrelation function K{t) is found by adding the three contributions ((67|) . (|78p . and (|82p . 

ii:(t) = X(<)(") + X(<)W + K{tY''\ (87) 

The classical probabilities P(C, B;t + t) and P(C, P; r) in Eq. (|67|) are the ballistic equivalent of the 'diffuson 
propagator' of the diagrammatic perturbation theory of universal conductance fluctuations; The time derivatives 
dP{B, A;TB)/dTB and dP{D ,C]t-e) / dr-^ are the ballistic counterparts of the Hikami box from diagrammatic per- 
turbation theory. A similar correspondence holds for the classical propagators in Eqs. ((78)) and (|82l) . although the 
roles of Hikami box and diffuson propagators are intertwined in these cases. Note that our expressions for K[t) do 
not contain the ballistic equivalent of the 'Hikami hexagon'. [A 'ballistic Hikami hexagon' appears, however, in the 
semiclassical calculation of tr [Si2S\2)^ which is not considered here.] 

Sofar we have not discussed a second set of three contributions to K{t), obtained by time-reversing the trajectories 
a2 and (32- In the language of diagrammatic perturbation theory this is the 'Cooperon contribution' to the conductance 
fluctuations. At zero magnetic field the formal expressions for the three Cooperon contribution to K{t) are identical 
to those of Eqs. ([S7)) . ([75)) . and (|82p . so that K{t) is increased by a factor two with respect to the expressions given 
above. At a finite magnetic field the classical propagation probabilities for the Cooperon contribution are multiplied 
with phase factors exp(2ie$/ft-c), where $ is the flux enclosed between ai and 02 between P and C [for iir(i)*^°^] or 
the flux enclosed by the periodic trajectory 7 [for K{t)'^^^ and K{t)'''^^]. 

The three contributions K{tY°'\ K{t)^''\ and Kit)'^'^^ have different Ehrenfest-time dependences. The first contri- 
bution, K{t)'^'^\ is finite and positive if te ^ and vanishes in the limit of large rg. The third contribution, K{t)'''^\ 
on the other hand, vanishes if te/td — > and grows upon increasing te- Finally, the TE-dependence of K{tY^') is 
non- monotonous: Kit)^^"^ equals its zero-TE limit if < = T^/n, n = 1, 2, 3, . . ., irrespective of the ratio te/td, whereas 
Kit)'^^'^ = for t = TE/{n — 1/2). These two remarkable properties follow from the fact that the integral in Eq. ([75]) 
is periodic in te with period t and symmetric around tAB.j — t{n — 1/2), j — 1,2, before taking the derivatives to 
tAB,i and tAB,2- Two examples of the te dependence of K{t) will be discussed in Sees. IV Cl and IVDl below. 



B. Encounters touching the lead opening 

The calculation of K{t) shown above was based on a semiclassical calculation of the covariance of the reflection 
coefficients off the two contacts. This method to calculate K{t) is technically simplest because one only needs to 
consider encounters that reside in the interior of the sample. The encounters that appear in the calculation of K{t)'^°''^ 
and K(t)^'^'> can touch the contacts, however, if a different expression is used to calculate K{t). The encounters that 
appear in the calculation of K{tY^^ never touch the contact because they are part of a periodic trajectory. In order 
to illustrate a semiclassical theory for K{t) with encounters that touch the lead opening, we calculate K{t) from the 
covariance of reflection and transmission coefficients, using 



g{e) = A^i -tr5ii(e)5ii(£)t, g{e') = tr 5i2(e')^i2(e')^ 



(88) 
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in Eq. ((60)) . In this case, all trajectories exit the sample through contact 1, so that the encounters may touch the lead 
openings upon exit. Upon entrance encounters still can not touch the lead openings because the trajectory pairs ai, 
Pi and a2, P2 enter through different contacts. For the three contributions to K{t) one then finds, in the absence of 
time-reversal symmetry 



dAdBdC P{A, 1)F(A, 2)—P{B, A; te) 



X I dTP{C,B;T + t)P{C,B,T) 



F(1,C;te)+ / dDP{l,Dy—P{D,C;TE) 

CITE 



i^W^'' = - J dAidA2dBidB2 dri2P(Ai,l)P(A2,2)P(l,Bi)P(l,B2) 



d 



d 



-Py{Ai, A2, Bi, B2; \t\,Ti2,tAB,l,tAB,2) 



tAB,l—iAB,2—'^¥: 



dtAB,l dtAB,2 

K{t)^^^ = - J dtBA J dAdBdCP{C,l)P{C,2)^P{A,C;TE-tBA)P^{A,B;\t\,tBA) 

X P{\,B-TE~tBA)+ ! dDP{\,Df^P{D,B-TE~tBA) ■ 



(89) 



(90) 



(91) 



One verifies that each of the three contributions separately equals the corresponding contributions ([67)) . ((78)) . and 
calculated using Eq. (p^ . 



C. Quantum dot 



As a first application of the general expressions for K{t) derived above, we consider the conductance fluctuations 
in a chaotic quantum dot. The probabilities P{A,j) = P{j,A) = Nj/N, j = 1,2, for all phase space points 
A in the interior of the quantum dot, while the probabilities P{A,B;t) of propagation in the dot are given by 
P{A, B; t) — il^^ exp(— i/TD), where f2 is the dot's phase space volume. With this, one finds 



if ^ illil2_g~2rE/TD-|t|/rD 



(92) 



For the calculation of K{tY ' we note that the probability P^{Ai^ A2, Bi, B2; \t\, T12, tAB.i, tAB.2) — ^ exp(— |i|/rD), 
independent of the propagation times tAB,i and tAB,2 between the phase space points Ai and Bi and A2 and B2. 
Since K{t)^^^ contains partial derivatives to tAB,i and tAB,2, one has 



K{tY^'^ = 0. 

Finally, the probability P^{A,B; \t\,tBA) = f^^^e^l'l/^n, so that 



(93) 



(c) _ NtN^ 



2 m2 i-te 
T) Jo 

2rE/TD^g-|t|/TD^ 



(1 



Adding Eqs. ([92]) and jM]), one findi 



23. 



Kit) = iXMe-l*IAo 



independent of the ratio te/tq. As a consequence, the variance of the conductance reads^^ 

vary = 



at zero temperature, and 



var g 



6Ttb iV4 



(94) 



(95) 



(96) 



(97) 
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if Ttd ^ ?i- In the presence of time-reversal symmetry, K{t) and var g are a factor two larger. The remarkable 
TE-insensitivity of var g was first reported in Ref. [l^ on the basis of numerical simulations of the quantum kicked 
rotator at zero temperature. 



D. Lorentz gas 

The second example we discuss is that of the quasi one-dimensional Lorentz gas. As discussed in Sec. IIIIDl for this 
example the phase space coordinates appearing in Eqs. (|67|) . ([78ll . and ([82|) can be replaced by the coordinate x along 
the sample. The probabilities P{A,B;t) and P{A,j) = P(j,A) are then obtained from a one-dimensional diffusion 
process, see Eqs. P5)) ~(|I5 |) . 

Equation ((78|) contains a double partial derivative of P-y . These derivatives act both on the exponent in Eq. ([43]) 
and on the theta function. We separate both actions by writing 

drPix, x; t) = 5{t)5{x - x) + dtP{x, x'; t), (98) 

where the first term comes from the derivative of a theta function 9{t), and the partial derivative dt in the second 
term acts on the probability P{x, x'; t) without the discontinuity caused by the appearance of the theta function 9{t), 

dtPix,x';t) = ^—0{t)y ^i^e-^'"'/^ sin^sin^. 

Ltu L L 

Taking the double partial derivatives according to Eq. (j98p and setting tAB,j equal to te, j = 1, 2, we then find 

dtAB,i9tAB.2P'riAl, ^2, Bi, B2; \t\,Ti2;TAB,l,TAB,2)\^^^ ^^^^ 
dridr2P^{Ai,A2,Bi,B2]\t\,Ti2]Ti,T2) 

Tl—T2—t 

- - ri2) + <5(ti2)]<5(Ai - A2)dr, {P(B2, A2; T2)P(Bl, B2; Ti - T2 + 0+)P(Ai, A2; \t\ - n) 

- P(i?i, Ai; ri)P(i?2, Pi; T2 - Ti + 0+)P(Ai, A2; \t\ - T2)} |^^^^^^^- . (99) 
where we abbreviated 

t = TE mod (100) 

Substituting Eqs. p5)) - (|45p we then find that the three contributions to the conductance autocorrelation function 
K{t) read 

1 p^P^ l*l/''"D / „ \ 

K{tf) = ^ycL(a2-^^')'(e-(l*l-l^*-l*")('''+^')/2^/.2,2 

4td L ' V '^D 

^ 2e-(l*l+|2t-|t|l)'^V2ro l^\t\-\2t-\t\\\^ I 

-- V C^.C^.Cp.CpJ(^2-^.2)(^2_^2)^-(|tH2tH*ll)(P^V)/2.o^^,^ 

TD I \ TD 

,\f,,2 r.2\(„2 J2\^f,,2 ,,2\(J2 , ,2 M ^-cr^ |2t- 1 1 1 1 /td -f ( 1^1 ~ 1^^ ~ 1^1 



/j,o- odd i^,p ^ 



+ /^2+p2+,2,,2 ( — j K . (103) 
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where 
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IF 2^ 




if /i^ 

else, 



- v odd, 

if /i + p 
else, 



odd 



if fl^ = 

else 



(104) 



In the presence of time-reversal symmetry, K{t) is a fac- 
tor two larger. 

Evaluating Eqs. pUT]) . pU^ . and piI5)) in the limit 
te ^ td one recovers the results 



-P^|t|/TD 



(105) 

(106) 

(107) 
disordered 



that are known from the theory of 

conductors .55,56 ,57 

We have not been able to evaluate Eqs. (jlOip . (|102p . 
and (|103p in closed form for finite te. The result of a 
numerical evaluation of the three contributions to K{t) 
for te/td = 0.2 and t-e/td = 2 is shown in Fig. [51 Note 
that K{t)^''^ is bounded from above by its zero-TE limit, 
and oscillates between zero for \t\ = T-E,/{n -f 1/2), n — 
0, 1, 2, . . . and the zero-TE limit for \t\ = T^/in + 1). 

The zero-temperature conductance variance is ob- 
tained by integration K{t) over t, cf. Eq. (pT|) . In the 
limit Te ^ td one thus finds the well-known result for 
the zero temperature conductance variance in a disor- 

5 

r 



dered quantum wire— '^'^'^ 



var 5 



1 

15' 



(108) 



In the presence of time-reversal symmetry, var g = 2/15. 
The TE-dependence of the zero temperature conductance 
variance is shown in Fig. [9] 

The conductance variance at high temperatures (not 
taking into account dephasing) is determined by the 
small-|<| asymptotics of K{t). The small-|i| asymptotics 
of K{t) is dominated by K^'^\ which diverges at small \t\, 



K{t) 



E 



/J, odd 



7' 



(109) 



Integrating over t, we find that the high-temperature 
asymptote of var g is given by 



var g = 



E 



1 



-2fj. te/to 



^ odd 

0{1/Ttu), 



(110) 
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FIG. 8: (Color online) The three contributions to the con- 
ductance autocorrelation function K{t) for te/td = 0.2 (top) 
and te/td ~ 2 (bottom). The dashed curves show K{t) for 
the limit te/td — > 0. 



where c « 0.0409 is a numerical coefficient. This is to be 
contrasted with the asymptotic temperature dependence 
of var g in the limit te/td — > 0, which reads 



1 



var 5 



0(1/Ttd) 



3/2 



(111) 



The small- |i I divergence of K{t)'^''^ comes about be- 
cause the return probability in a quasi one-dimensional 
Lorentz gas diverges cx t"-^/^ a,t short times (but times 
longer than the mean free time). This enhances the 
contribution to K(t) from trajectories that involve pe- 
riodic orbits compared to the contribution from trajecto- 
ries that do not. The second contribution that involves 
periodic orbits, K{t)'^^\ has an additional time integra- 
tion, since it requires that two encounters are placed on 
the same periodic orbit. This limits the available phase 
space at small \t\, so that K{tY'''> oc It/Tnl^^^ for small 
\t\. 
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FIG. 9: Conductance variance versus Ehrenfest time for a 
quasi one-dimensional Lorentz gas. The figure shows var g 
in the absence of time-reversal symmetry. In the presence of 
time-reversal symmetry, var p is a factor two larger. 

VI. DENSITY OF STATES 

It is instructive to compare our calculation of the meso- 
scopic fluctuations of the conductance of an open quan- 
tum system to a calculation of mesoscopic fluctuations of 
the density of states. 

Fluctuations of the density of states are described by 
the spectral form factor K^{t). The definition of the 
spectral form factor is similar to that of the conductance 
autocorrelation function in Eqs. ((59|) and ([60| . 

K^{t) = — f dLuK„{hLj)e''^* (112) 
2vr J 

K^e-e') = {^{e)v{e'))-{v{e))\ (113) 

The trajectory-based semiclassical evaluation of the spec- 
tral form factor in a closed quantum system commonly 
starts from the Gutzwiller trace formulaj^i^ 

K,{t) (X AtA'^^e^^^^^~^fyH[t (r„ + r^)/2]\ , 

\a,/3 / 

(114) 

where a and fi are periodic orbits of duration Tq, and T/j, 
respectively. (See Ref. [6l| for a field-theoretic calculation 
of K^{t).) Here A'^' and S'^^ are the stability amplitude 
and classical action of the periodic trajectories in the 
closed system, respectively. Their precise definitions dif- 
fer slightly from the definitions of the stability amplitude 
and classial action used in the semiclassical expression of 
the scattering matrix of an open quantum system, see 
Sec. nil hence the superscript "cl" . The leading contri- 
bution to Ki,{t) in the semiclassical limit Ti ^ comes 
from the diagonal terms a = (3 (up to time-reversal, if 
time-reversal symmetry is present) The multiplicity of 



each periodic orbit is proportional to its duration. In the 
language employed in this article, one then findaii 

^"^^'^ ^ (2^ / '^^l^l^-rC^; \t\). (115) 
where the integral extends over all phase space points 
A and P.y{A; \t\) is the probability that A is on a peri- 
odic trajectory with period \t\. In the presence of time- 
reversal symmetry K^{t) is a factor two larger. Setting 
P~f = ri~^ Eq. (|115p reproduces the form factor of ran- 
dom matrix theory;®*' Taking P-y{A; \t\) to be the return 
probability in a diffusive medium one recovers the spec- 
tral form factor of a disordered metal. ''^'^^ Since no tra- 
jectories with small-angle encounters are involved in the 
semiclassical calculation of K^{t), the leading contribu- 
tion to K^{t) has no Ehrenfest-time dependence. (This 
is different for quantum corrections to that are of 
higher order in Ti^iL^ Such quantum corrections are not 
considered here.) 

In an open quantum system, the density of states is 
formally expressed in terms of the Wigner-Smith time- 
delay matrix, 

= -il^^'-^'^^^^^^ + ^^^kio- (116) 

Nevertheless, the fluctuations of the density of states can 
still be calculated from the Gutzwiller trace formula, in 
which case the trajectory sums are naturally cut off for 
Ta, Tp ^ Tu- The result (|115p for the leading contribution 
to K^{t) also applies to open conductors. 

If, on the other hand, the fluctuations of the density 
of states are calculated using Eq. ()116p . one arrives at 
a summation over sets of four classical trajectories, not 
pairs. A priori, this summation is very similar to the sum- 
mation over sets of four trajectories needed to calculate 
conductance fluctuations. Yet, the final result must equal 
the simple expression (jllSp obtained from the Gutzwiller 
trace formula. We now show how this simplification oc- 
curs. 

Using the semiclassical expression for the scattering 
matrix, wc write 1^(6) as 

(117) 

where the trajectories a and /3 enter (exit) the conduc- 
tor through the same contact and with the same stable 
(unstable) phase space coordinate Sq, {u'^). 

Before we calculate the fluctuations of the density of 
states, it is instructive to first verify that there is no 
interference correction to the average density of states. 
This can be seen explicitly from the results of Sec. IIV| 
which give 
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Sv{e) = TTT f dtdt'dt" I dBdC^{t' +t" + 2TE + t)P{B-t')P{B-t")P{C,B-T^)P{C,C-t) 
2nn J J dTE 

+ ^JdtJ dC{2TE + t)P(C; m)P(C, C; i), (118) 

where P{B\ t) — P(l, B; t) + P(2, B] t) is the probabiUty to reach any contact from the phase space point _B in a time 
t. Using the identity 

j dtP{B-t) = 1, (119) 

and Eq. ([55]) this is rewritten as 

5v{e) ^ dtdt' J dBdC (^^'^ + l) P(B; t')P{C, B; te)P(C, C; t) 

= ^J'^idt'^ J dBdCt'P(B;t')PiC,B;TE)P(C,C;t) 

- 0, (120) 

in agreement with our expectation that there be no quantum corrections to the average density of states. 

For the calculation of the spectral form factor K^{t) one finds three contributing configurations of classical tra- 
jectories. These are the same as the three configurations of classical trajectories that contribute to the conductance 
autocorrelation function K{t), see Fig. [6^-c. The contributions from the trajectory configurations of Figs. [6^ and c 
vanish, however, by arguments similar to those used in Eqs. pisp and (jl20p above. Hence, we only need to consider 
the trajectories of the type shown in Fig. [Bb, for which we find 

1*1 



Ky{t) = J^^^ J dAidA2dBidB2 J dt'^dt'^dt'ldt'^ dns 



^ ^ -{t[+t'( + tAB,l + \t\){t'2+t'^ + tAB,2) 



dtABS dtAB,2 

X P(3i;i;)P(A2;i;,)P(Bi;t'/)^(S2;i2)^7(^i.^2,Si,B2;|t|,Ti2,iAS,i,tAS,2) 



■ (121) 

tAB,l—tAB,1—T¥. 



Either P{Ai\t\) or P{Bi;t'{) will be integrated over the 
time argument and, using Eq. (|119p . drop out of the ex- 
pression for Ky{t). Since P^ no longer depends on tAB,i 
once integrated over Ai or Bi, the partial derivative to 
tAB,i must act on the prefactors t'^ -\- 1'( + Iab.i + \t\, not 
on P^. Similarly, the partial derivative to tAB,2 must act 
on the prefactor ^2 + *2 + "^e- Again using Eq. pi9p . we 
remove the remaining two of the four factors P(i?i;i"), 
P(Ai;t[), P(A2;t2), and P(i?2;t2) from the expression 
for K^{t). With this, we find precisely the result that fol- 
lows from the Gutzwiller trace formula, Eq. (|115p above. 

It is interesting to point out that, while there was 
a simple reason why K^{t) does not depend on the 
Ehrenfest time — in the Gutzwiller trace formula Ku{t) 
is calculated from trajectories without small-angle in- 
tersections — , the corresponding contribution K(t)^'''> 
to the conductance fluctuations generically has a te- 
dependence, as is illustrated, e.g., in Fig. [51 Also, we 
note that the second contribution to the conductance 
fluctuations that survives in the limit of large Ehrenfest 



times, K^''\t), has no counterpart in the theory of spec- 
tral fluctuations, even at finite Ehrenfest time. This is 
particularly remarkable for the case of a chaotic quan- 
tum dot, where K''^\t) = 0, so that the density of states 
fiuctuations and the conductance fluctuations at large te 
are due to two different conflgurations of classical trajec- 
tories! 



VII. CONCLUSION 

In this article we derived expressions for the Ehrenfest- 
time dependence of three signatures of quantum trans- 
port: The Fano factor for the shot noise power, the weak 
localization correction to the conductance, and the con- 
ductance autocorrelation function. Our derivation is for 
arbitrary ballistic conductors with chaotic classical dy- 
namics and makes essential use of the inequalities Xp <^ I 
and r <C td, where Xp is the Fermi wavelength, I = vpT 
is the mean free path, and td is the typical dwell time 
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for electrons that contribute to transport. We did not 
take into account the effect of a time-dependent bias or 
the effect of dephasing from a time-dependent potential 
inside the sample. Our results agree with known results 
for the weak localization correction and Fano factor in 
arbitrary ballistic conductors,— and for conductance 
fluctuations in ballistic quantum dotsi^i^ 

The Ehrenfest-time dependence of the conductance 
fluctuations is remarkable, not only because the con- 
ductance fluctuations remain finite in the limit of large 
Ehrenfest times, but also because there are significant 
differences between the conductance fluctuations in bal- 
listic conductors at large te and in disordered conductors. 
These qualitative difference come from of the existence 
of a configuration of classical trajectories that contribute 
to the conductance fluctuations in ballistic conductors 
at flnite te, but not in disordered conductors. This 
class of trajectories is shown schematically in Fig. [6t. 
One important difference is the temperature dependence 
of the conductance fluctuations: In disordered quantum 
wires, var g is proportional to for large temperatures, 
Ttd ^ h, where td is the typical dwell time for elec- 
trons travelling between source and drain contacts. If the 
point-like impurities in the quantum wire are replaced by 
macroscopic discs, asymptotically var 5 cx This 
asymptotic temperature dependence should be visible for 
temperatures above Titd/ min(rE,TD)^. 

Our final results arc formulated in terms of classical 
propagators that represent the possibly correlated prop- 
agation of quantum-mechanical transition amplitudes 
along up to six classical trajectories. It is relatively 
straightforward to include a time-dependent bias or the 
effect of slow variations of the potential and/or mag- 
netic field inside the sample in the theoretical framework 
presented here. A time-dependent bias requires one to 
consider products of the sample's scattering matrix S at 
different energies. ^'^ Since small and spatially slow varia- 
tions of an internal potential or a magnetic field do not 



alter classical trajectories, their effect can be included by 
a suitable modification of the classical action Sa for each 
trajectory a.— 1^ Time-dependent internal potentials can 
be treated in a time-resolved semiclassical approach^ In 
all cases, the only change to the theory as it is presented 
here is that one has to weigh the classical propagators 
with an additional phase factor representing the accu- 
mulated action differences. The TE-dependence of the 
conductance fluctuations in a ballistic quantum dot in 
the presence of a magnetic field were considered in Ref. 
[23. 

The situation is more complicated if one considers 
a potential internal to the sample that varies on sub- 
macroscopic length scales. Such a potential may arise 
as a result of electron-electron interactions inside the 
sample,— or it may represent residual point-like impu- 
rities. The complications arise because such a potential 
may alter the electron's trajectory. (Potentials with spa- 
tial variations on macroscopic scale cannot change the 
electron's momentum by more than Ti/l. Such small mo- 
mentum changes have no consequences for times ^ te.) 
Petitjean et al. argued that a trajectory-based semiclas- 
sical approach still can be used as long as the potential 
is smooth on the spatial scale ~ [IXpY^'^. How to in- 
clude potentials with a faster spatial dependence into the 
present formalism remains an open question. 
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